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Abstract 

We present results for total cross sections, single and double differential distributions 
and correlations between pairs of outgoing particles in the reactions p + p ^ + 7 and 
p + p ^ + 7 + jet at -sfS = 1.8 TeV. Order as QCD corrections and leading logarithm 
photon bremsstrahlung contributions are included in the MS mass factorization scheme for 
three experimental scenarios: 1) 2-body inclusive production of and 7, 2) exclusive 
production of , 7 and 1 jet and 3) exclusive production of and 7 with jet. The 
latest CTEQ parton distribution functions, which fit the newly released HERA data, are 
used in our analysis. The dependence of our results on the mass factorization scale is used to 
place error bars on our predictions for the single differential distributions and correlations. 
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I INTRODUCTION. 

In the first part of this study [I]] we presented the general formahsm for the computation of 
exclusive cross sections in p + p —>■ + 7 and p + p —>■ + 7 + jet. We examined the 
order as and photon bremsstrahlung contributions to these cross sections in the MS mass 
factorization scheme. The analytical expressions obtained in that work were then combined 
with a Monte Carlo integration routine, namely Vegas [3, to produce numerical results 
for three different experimental scenarios: 1) 2-body inclusive production of and 7, 2) 
exclusive production of W^, 7 and 1 jet and 3) exclusive production of and 7 with jet. 
The resulting programs have been run for ^/S = 1.8 TeV in the proton-antiproton center 
of mass system (Fermilab Tevatron). We include 4 massless flavours u,d,s and c in our 
partonic hard scattering processes and set cos 9c = 0.95 where 9c is the Cabbibo angle. We 
use the one loop strong running coupling constant with Aqcd = 0.139 GeV. The W boson 
is assumed to be an on shell physical particle with mass = 80.2 GeV. 

We use the latest available CTEQ parton distribution functions in the MS mass 
factorization scheme (set CTEQ2M). These parton distribution functions are in agreement 
with the latest CCFR NLO analysis of the strange quark density and flt the newly released 
HERA data (for details see 0.) 

The three experimental scenarios considered here are deflned as follows (see also the dis- 
cussion in m): 



1) 2 BODY INCLUSIVE PRODUCTION OF AND 7. 



In this scenario ("2 body inclusive scenario") the reaction p+p W'^ +'~f+X with X = 
or 1 jet is considered. In other words, one detects the outgoing boson and the photon 
but does not tag the outgoing jet. We therefore use the following kinematic restrictions on 
the outgoing particles, in the proton-antiproton center of mass frame: 

cos 6^^ I , I cos 6^^^/ 1 < cos(0.3 rad) 

10.0 GeV 
0.2 

< 0.2) =^ (S0et,7) < 0.2) 
(%t,w^<0.2) =^ (s(,-,t,M^) < 0.2) , (1.1) 

where 9i is the angle between the incoming proton axis and the axis of the outgoing particle 
i and Pti is the transverse momentum of particle i. Rij is the cone size between a pair of 



, COS 9w 


< 


Pt^, Pt\Y 


> 




> 



outgoing particles: Rij = \J{ri*.jY + {(t)i,jY with rj* the pseudorapidity and the azimuthal 
angle; S(^jet,w) = Ejet/Ew is the "shadowing ratio" between the untagged jet and the 
boson. The third condition in ( |1.1| ) removes events where the boson is too close to the 
photon. The last two conditions remove events where the jet which is too close to the 
boson or photon is at the same time of comparable energy so that it would "shadow" one of 
the two tagged particles making it undetectable. This might be the case for bremsstrahlung 
contributions with a small photon momentum fraction. 

Several artiflcial parameters, namely xo,yo and vq were introduced in |l| in order to 
control the numerical cancellation of singularities. Results for physical quantities like total 
cross sections and differential distributions and correlations do not depend on the choice of 
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these parameters, however, the requirement of stabihty and small numerical errors in the 
Monte Carlo runs lead us to the choice xq = (1 + p(s))/2, i/q = 1.0 for the partial cross 
sections in the qq channel and yo = Vo = 0.05 for the partial cross sections in the qg and gq 
channels. For the cut parameters defined in section V of |1|] we choose the following values: 
Ax = 10~^ and Ay = A^ = 10^^. As discussed in the referred section these cuts introduce 
errors in the numerical results, however, we have checked that the above choice minimizes 
these errors while keeping the Monte Carlo program stable. Our tables of results for total 
cross sections include entries for the 

lowest order estimates of these errors and confirm that they are negligible. Thus the 
expressions for the errors are neglected in all the distributions and correlations. 



2) Exclusive production of 7 and 1 jet. 



In this scenario ("1 jet scenario") one detects three outgoing particles, namely the W^, 
7 and 1 jet. This scenario is thus defined by the following kinematic conditions: 



cos 9w , cos 9 jet 


< 


cos(0.3 rad) 


Pt^, Ptwi Ptjet 


> 


10.0 GeV 




> 


0.2 


Rjetj-y 


> 


0.2 


Rjet,W 


> 


0.2. 



The above conditions will automatically remove bremsstrahlung contributions which tend to 
produce jets parallel to the outgoing photon. In the Monte Carlo runs we choose the values 
for the artificial parameters of this scenario in all channels as: xq = 1, yo = vo = 0- This 
choice and the above cuts guarantee that no unphysical dependence is introduced in the cal- 
culation of physical quantities related to the jet. Moreover, the cuts A^., Ay, A^, introduced 
in the first scenario are not necessary here and thus the associated errors are zero in this case. 



3) Exclusive production of and 7 with jet. 



In this scenario ("0 jet scenario") we select events where the boson and the photon 
are detected and no outgoing jet is detected. This will include 2 to 3 body events where the 
outgoing jet is not detected because it has either a small angle with respect to the beam 
or a small transverse momentum or it is "shadowed" by the W'^ boson or the photon. The 
following kinematic conditions define this scenario: 

I cos 6^^ I , I cos I < cos(0.3 rad) 



, COS 9w 


< 


Pt^, Ptw 


> 


Rw,'^ 


> 



10.0 GeV 
0.2 

{Rjet,^ < 0.2) =^ {S(jet,^) < 0.2) 
{Rjet,W < 0.2) =^ {S(jet,W) < 0.2) 

(I cos9jet\ > cos(0.3 rad)) or {Ptjet < 10.0 GeV) . (1.3) 

As already discussed in [|I| the results for this scenario can be obtained by subtracting the 
results for the 1 jet scenario from the corresponding results for the 2 body inclusive scenario. 
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In section II we present results for total cross sections. Section III is devoted to the 
discussion of single differential distributions and correlations. In section IV we present and 
discuss results for double differential cross sections and correlations. In section V we end 
our study with some concluding remarks. 
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II TOTAL CROSS SECTIONS. 
II.A Brief Review of Relevant Formulae. 

We write the total hadronic cross section as: 
Jo Jo 

+ (p ^ P,Ti ^ T2,Pl ^ P2, Pi P2)} 



(2.1) 



with 



Born P P 

^qq = ^qq + ^<?<? {SV) + ^la, + lb + 1 ,4, + CTgg (finite) + ^99 {Brems) + Cfqq^error) 

— rr^'^ -I- n-^'™' -I- rr^'^^ -I- rr^^^'^"^ ?rP rr -I- 

"^"JS ^qg, finite "^gg ^qg, finite "^gg "^qg ^qg(Brems) "T C'qg{error) 

_ P,/ I _./,co« I _|_ -.//,coi I =-P I ^ _ _l_ n- - 

"^OT ~ "^gq, finite "^gg '^gq, finite "^gij "-"gg ^gq(Brems) ~T (Jgq{error) 



(2.2) 



and 



C/b = <^Ib,l + C'"7'fe,2 + 0"/6,3 

(^qq {finite) = ^f,l,lA + Cr/,l,2,a + Cr/,l,3,a + Cr/,l,l,fc + ^f,l,2,b + 0-/,l,3,fe 

+ Cr/,2,l,a + C'"/,2,2,a + f,2,3,a + '^/,2,l,fe + '^/,2,2,fe + f,2,3,b + 0"/,3 

P,/ — rr^ -L rr^ J- rr^ 

^qg, finite " '^gg,/,l '^9g,/,2 '^gg,/,3 

Oqg, finite " ^gg,/,l + '^55,/,2 + ^gg,/,3 

P,/ _ 7 . I . I 

^gq, finite " ^^gg,/,! + ^gqj,2 + '^gg,/,3 

^gqjinite " ^gg,/,! + '^(;g,/,2 + Crgg,/,3 • l^-^J 



Each term in ( p.2|) and ( |2.3| ) has been explained in section V of [|1|. Note that the 
replacements in parenthesis in ( ^.1[ ) apply to all explicit as well as implicit variables in the 
above expressions. The last replacement in the parenthesis in ( p.l| ) will act on physical 
variables that go into experimental cut functions and histograms (see section IV. D of [|l|.) 

For the bremsstrahlung contributions we use here the so called leading logarithm approx- 
imation, which has been previously used by other authors In this approximation 
we write the photon fragmentation functions as: 

M2 



fiq(q)i.X,M) 




A2 



e2(-) (2.21 - 1.28x + 1.29x2) x'^-^^^ 
1 - 1.631n(l -x) 

+0.002(1 - x)2-0x-2-54l 
f^g{x,M) = —In ( 1 0.0243(1 -x)^■°^x-^■^^ (2.4) 




where eg(g) denotes the charge of the outgoing quark (antiquark) q (g) in units of e, M 
is the mass factorization scale and the running electromagnetic fine structure constant is 
a = e2(/i2)/47r . 
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II.B Total cross section for the 2 body inclusive reaction p- 



All terms in expressions (|2.2| ) and ( p. 3D will produce nonvanisliing contributions in this 
scenario. In table ^ we show the results for the qq channel cross sections in pb at the three 
scales r = 0.5, 1.0 and 2.0, where r = M/Mw = ii/My/ and M, are the mass factorization 
and renormalization scales, respectively. The second entry of table shows the value used 
for as at these scales. In all our runs we use a = 1/137.036 for the fine structure constant 
and Gp = 1.16639 ■ lO^^GeV"^ for the Fermi coupling constant. 

The corresponding results for the qg channel contributions to the 2-body inclusive hadronic 
reaction p + p —>■ + 7 + X are summarized in table |^. Finally, the corresponding results 
for the gq channel are summarized in table |^. 

We thus obtain total cross sections for the 2-body inclusive reaction p + p ^ + 7 + X 
of 4.70,4.28 and 4.04 pb for r = 0.5,1.0 and 2.0 respectively. The total cross section at 
order as is thus a decreasing function of the scale r, confirming the results obtained in a 
previous study |^. The strong dependence of the cross section on the scale r means that 
order a| corrections are large and they must be included for accurate predictions in this 
scenario. This was also pointed out in |^ when studying the single particle photon inclusive 
cross section. 

If we compare these results with the ones previously reported for the single particle 
inclusive total cross section in |^ we note a considerable reduction of the total cross section 
in the present calculation. The smaller values can be attributed to larger cut in the photon 
angle and the effect of experimental cuts on the boson which was treated previously in 
an inclusive fashion. 

We also note that in the approximation we are using here the bremsstrahlung contri- 
butions are negligible in the qq channel and that their contribution in the other channels 
amounts to not more than 2% of the total cross section. This would make big deviations 
from the leading logarithm bremsstrahlung approximation easy to study. 

The next to last rows in tables and |^ show that the errors introduced by the A^., Ay 
and A^ cuts (see section V of |jl|],) are negligible and thus our prescription for splitting the 
x,y and v integrals when adding up the histograms is consistent. 



II. C Total cross sections of p+p W^+^+jet and p+p 

With the choice of parameters Xo,?/o and Vq explained in section I the only terms in ( p.2| ) 
contributing to the reaction p + p + 7 + jet are: 

= ^^OT-j.i + ^Sja + ^M- (2-5) 

The results for the contributions of all these terms to the hadronic cross section are shown in 
table ^ for the three scales r = 0.5, 1.0 and 2.0. At these scales we thus obtain for the total 
cross sections of the 1 jet reaction values of 2.15, 1.79 and 1.51 pb respectively. As in the 2 
body inclusive case the variation is large in the 1 jet scenario so we need to include higher 
order QCD corrections to make more accurate predictions. Unfortunately these higher order 
corrections are not available. 
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r 


0.50 


1.0 


2.0 


a.s 


0.145 


0.129 


0.116 


qq 


2.66 


2.62 


2.58 


' .SI 


0.31 


0.51 


0.67 




n 09 
u.uz 


u.uz 


n 09 
u.uz 


C/a,2 


u.uo 


-u.uo 


-U. i i 


C/a,3 




1 8 9^; 
-io.zo 


- i ( .yo 




-18.39 


-18.28 


-18.07 




0.02 


0.02 


0.01 




0.04 


-0.04 


-0.10 




-18.13 


-18.23 


-17.91 




-18.38 


-18.25 


-18.03 




37.20 


36.83 


36.29 


(^f,l,l,a 


0.07 


0.06 


0.05 




1.41 


1.20 


1.04 


0-f,l,3,a 


-1.30 


-1.12 


-0.97 




0.08 


0.07 


0.06 




1.20 


1.03 


0.90 




-1.14 


-0.98 


-0.85 




1.12 


0.98 


0.87 


'^/,2,2,a 


28.97 


26.45 


23.44 




-29.03 


-26.51 


-23.50 




1.12 


0.99 


0.87 


<^/,2,2,6 


29.95 


26.31 


23.33 


<^/,2,3,6 


-30.01 


-26.37 


-23.38 




-2.30 


-2.02 


-1.79 


p 

^gq(finite) 


0.16 


0.10 


0.07 


^ qq{Breras) 


0.003 


0.003 


0.003 


^qq{error) 


0(10-4) 


0(10-4) 


0(10-4) 


(Jqq 


3.61 


3.52 


3.52 



Table 1: Results for the partial hadronic cross sections (in pb) for the 2-body inclusive 
reaction p -|- p — > -|- 7 -|- X in the qq channel at =1.8 TeV. 
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r 


0.50 


1.0 


2.0 


as 


0.145 


0.129 


0.116 


^aa f 1 


0.38 


0.31 


0.26 




0.78 


0.66 


0.56 




-0.80 


-0.67 


-0.57 


^qg, finite 


0.38 


0.31 


0.26 


^l,col 
^qf) 


-0.22 


-0.25 


-0.26 




0.12 


0.10 


0.09 




0.15 


0.13 


0.11 




-0.08 


-0.07 


-0.06 


i> 1 1 

r^i 1 
^qq, finite 


0.19 


0.16 


0.13 


^11 , col 


-0.04 


-0.04 


-0.04 




0.15 


0.12 


0.10 


(^qg(Brems) 


0.02 


0.02 


0.02 


(^qg{error) 


0(10-9) 


0(10-9) 


0(10-9) 


(^q9 


0.48 


0.33 


0.21 



Table 2: Results for the partial hadronic cross sections (in pb) for the 2-body inclusive 
reaction p + p ^ -|- 7 -|- X in the qg channel at y/S =1.8 TeV. 



r 


0.50 


1.0 


2.0 


as 


0.145 


0.129 


0.116 




0.14 


0.12 


0.10 




0.60 


0.50 


0.42 




-0.61 


-0.51 


-0.43 


qq, finite 


0.13 


0.11 


0.09 


I ,col 

(^gq 


-0.18 


-0.20 


-0.20 


S" 


0.38 


0.31 


0.26 


^qq,f.2 


0.48 


0.40 


0.34 


^qqj,i 


-0.27 


-0.23 


-0.20 


P,ll 

qq, finite 


0.59 


0.49 


0.41 


ll,col 

^gq 


-0.13 


-0.13 


-0.13 


~P 

^qq- 


0.11 


0.09 


0.07 


^ gq{Brems) 


0.08 


0.08 


0.07 


'-^ gq(error) 


o(io-«) 


o(io-«) 


0(10-*^) 


^gq 


0.61 


0.43 


0.31 



Table 3: Results for the partial hadronic cross sections (in pb) for the 2-body inclusive 
reaction p -h p — > -t- 7 -t- X in the gq channel at =1.8 TeV. 
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r 


0.50 


1.0 


2.0 




0.145 


0.129 


0.116 








n '?Q 


(7 f ^ 1 h 


0.52 


0.44 


0.38 




1.05 


0.89 


0.77 




u.zu 


n 91 


n 1 7 


Jl 


14 


1 1 


OQ 




0.12 


0.10 


0.08 


'III 


0.52 


0.12 


0.35 




0.07 


0.06 


0.05 




0.41 


0.33 


0.27 




0.10 


0.08 


0.06 




0.58 


0.47 


0.39 


^ total 


2.15 


1.79 


1.51 



Table 4: Results for the partial and total hadronic cross sections (in pb) for the reaction 
p + p ^ + 7 + jei at = 1.8 TeV. 

Subtracting the above numbers from the ones for the 2 body inclusive reaction p + p — 
W^"*" + 7 + X we obtain for the jet reaction p+p — > W^"*" + 7 the values of 2.55, 2.50 and 2.53 
pb at the scales r = 0.5, 1.0 and 2.0 respectively. In this scenario we do not need to include 
order q;| QCD corrections to make our predictions more reliable, since they are already very 
stable. 
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III SINGLE DIFFERENTIAL CROSS SECTIONS AND 
CORRELATIONS. 

We now turn to the analysis of the single differential distributions and correlations. Our 
results are shown in figures 1 to 25. The error bars represent the theoretical uncertainty 
associated with the dependence of our results on the mass factorization and renormalization 
scales. They have been obtained evaluating the distributions at the two scales r = 0.5 and 
2.0. Note that the central values and the upper and lower limits still contain a numerical error 
introduced via the Monte Carlo. This error is negligible in our 1 jet results but it is larger 
in our 2 body inclusive and jet results. Regions with big error bars in our plots may be 
interpreted as regions where perturbation theory at order as is not reliable and thus higher 
order QCD corrections would be needed to make more accurate theoretical predictions. 

In general, the single differential distributions and correlations have very little dependence 
on the scale r in the scenario where the and photon are produced with outgoing jet. 
If we include contributions from 1 jet processes we increase the statistics by more than 
60%, but a nonnegligible theoretical uncertainty is introduced by the scale dependence of 
the 1 jet processes. To make theoretical predictions for single distributions and correlations 
more reliable in the 2 body inclusive scenario, we would thus need to include the order q;| 
corrections. 

Figures 1 to 6 show results for single differential distributions for the boson and 
figures 7 to 11 show the corresponding photon distributions. These plots demonstrate that 
the W boson is mainly concentrated in the regions of small energy and transverse momentum 
in all three scenarios. In the 1 jet scenario peaks at around 85 GeV (figure 1) while 
Pty/ peaks at around 30 GeV (figure 2). In the other two scenarios the single differential 
distributions of Ew and Ptw are monotonically decreasing functions except for the physical 
and experimental lower cuts. 

Figures 3, 4 and 6 show clearly that there are dips in the polar angle and pseudorapidity 
distributions of the W boson at around 9w = 7r/2 in the jet and 2 body inclusive scenarios, 
with smooth dips at around the same point in the 1 jet scenario. This is very different to 
what we observe for the photon polar angle and rapidity distributions, which present dips 
at around 9w — 27r/3 in the jet scenario only (figures 9, 10 and 11). 

The jet scenario would so far be the best for studying deviations from the Standard 
Model by looking at photon distributions since on the one hand the theoretical uncertainties 
related to dependences on the scale r are small and on the other hand the 2 to 3 body 
contributions are suppressed so that the partonic zeros still show up as wide dips in the 
angular and rapidity distributions of the photon (see figures 9, 10 and 11). The 9-y distribution 
maintains its features when 1 jet events are added, as seen in figure 9, however, the dip in 
rapidity is smeared out, as seen in figure 11. If photon bremsstrahlung contributions are not 
well accounted for by the leading logarithm approximation, these may disturb the dips in 
the photon angle and the photon rapidity and all other photon single and double differential 
distributions and correlations in the jet scenario. Therefore, to isolate this effect it may 
be important to look also at physical quantities in the 1 jet scenario, where 

there is no contribution from photon bremsstrahlung at present order in perturbation 
theory. 

Figures 12 to 15 show single differential distributions of the outgoing jet in the 1 jet 
scenario. We see that the jet is concentrated in the low energy, low Pt and central rapidity 
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regions. We also note that the scale dependence increases in these regions. The plot in figure 
14 shows a quite uniform polar angle distribution for the jet. 

Figures 16 to 25 summarize our results for single correlations between outgoing pairs 
of particles. The same comments about the scale dependences in the different scenarios 
observed in figures 1 to 15 are valid here. 

The Rw,-y cone size correlation in figure 16 presents a sharp peak followed by a sharp 
dip that falls below zero near Rw,'y = vr in the jet and 2 body inclusive scenarios. A 
similar feature was observed by Mangano et al. in Rg g for heavy quark pair correlations 
with fixed quark rapidity. The corresponding correlation in the 1 jet process presents no 
anomalous behavior, suggesting that the observed dips in the other two scenarios are the 
effect of the subtraction pieces introduced when the factorization theorem is implemented to 
cancel coUinear singularities associated with untagged jets produced parallel to the beams. 

The Rwjet cone size correlation in figure 21 has a very symmetric shape with a peak 
centered at vr. On the other hand, Rjet,'y in figure 21 is much smoother between 0.2 and 
TT, falling rapidly to zero outside this range. This means that the jet and the photon are 
approximately uncorrelated in cone size, however the jet and the W boson are highly corre- 
lated. Figure 25 confirms that the latter pair of particles are mainly concentrated on a plane 
that contains the beams, i.e. the distribution of the azimuthal angle difference 4>w,jet peaks 
at TT. However, (f)jet,-y is approximately flat over the whole range. The observed bump near 
(pjet,-y = Rjet,-y = 0.2 is causcd by the required experimental cut on the cone size. 

Pair mass correlations are of interest when studying deviations from the Standard Model, 
as pointed out in [|]. The pair mass correlations Mwj = {Qi + PjY with j = 7, jet in 
flgures 17 and 22 both peak at around 100 GeV, however the correlation Mjet^-y has a plateau 
between 20 and 40 GeV slowly falling to zero above 40 GeV. 

Figures 18 and 23 show the distributions of the angles aij between pairs of particles 
As expected aw,-y peaks near a = vr in the jet scenario, since here the main contributions 
come from back to back 2 to 2 body partonic reactions. The smearing is introduced when 
boosting the partonic system into the hadronic system. This is not the case for the corre- 
sponding correlation in the 1 jet scenario, which shows that the W boson and the photon 
are uncorrelated in the angle difference. The angular correlations between the jet and the 
W or the photon are also quite smooth. 

Since transverse momenta are preserved under Lorentz boosts along the beamline, the 
distribution of the azimuthal angle difference 0vy,7 only makes sense in the 1 jet scenario, 
which is shown in flgure 19. It also peaks at around 0vk,7 = meaning that the W and the 
photon are mainly concentrated on the plane that contains the beams. This, along with the 
information contained in flgure 25, shows that the beams and the three particles produced 
in the 1 jet scenario are approximately coplanar. 

Figure 20 shows the distributions of the pseudorapidity difference between the W and 
the photon in the three scenarios. In all cases the events are concentrated along the negative 
region, which means that the photon is produced with an angle with respect to the incoming 
antiproton beam that is usually smaller than the corresponding angle of the W boson. We 
note the presence of two dips in the jet scenario: one around f]\Y^ = —1 and the other 
around 77^^ = 1.3. These dips are smeared into plateaus in the inclusive scenario by the 
effect of 1 jet processes. In flgure 24 we see a very different behavior for the pseudorapidity 
differences of jet and W or photon: ?7^ ,,et is symmetric around while f]jet,-y receives bigger 
contributions on the negative axis but both correlations peak at 0. 
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IV DOUBLE DIFFERENTIAL CROSS SECTIONS AND 
CORRELATIONS. 

Double differential cross sections and correlations help us visualize the qualitative features of 
the spatial distributions of particles hitting an experimental detector. We present in figures 
26. a to 45 results obtained by averaging Monte Carlo runs at r = 0.5 and r = 2.0. Figures 
26. a to 27.C correspond to double differential distributions of the W while figures 28. a to 
29. c show the corresponding double differential distributions of the photon. Figures 30. a to 
34. c contain the plots of the correlations between the W boson and the photon and figures 
35 to 45 the various correlations between the jet and either the W or the photon. 

First we look at the double differential distributions of Pt vs 9 for the W boson (figures 
26. a, b and c). We conclude that the dip in 9w, already observed in the single differential 
distribution (figure 3), increases in relative depth as Pt^ decreases. In the 1 jet scenario the 
double differential distribution clearly peaks along the line Ptyy = 30 GcV for all values of 
9w (figure 26. b). In the other two scenarios the peaks move down to small values of Ptw, 
but in figures 26. a and 26. c they are not clearly seen due to the low density of bins that have 
been used on each axis; these peaks show up in the single differential distributions in figure 
2, however, where the density of bins is four times bigger. Similar conclusions are obtained 
when looking at Pt vs pscudorapidity for the W in figures 27.a,b and c. 

The corresponding double differential distributions for the photon have been plotted in 
figures 28. a to 29. c. The effect of the 1 jet processes on the dips in angle and rapidity 
which reflect the partonic zeros is clearly visible in these flgures. We note that in the 
jet scenario the dips increase their relative depth as Pt^ decreases. This suggests that the 
effects of deviations from the Standard Model on the single differential distributions in the 
polar angle and rapidity of photon would be enhanced if Pt^ is constrained to take values 
below certain bound, say 25 GeV. 

From the results in the previous subsection we ah^cady know that the events are mainly 
concentrated in the regions of small E]y and small E^^. Moreover, figures 30. a, b and c show 
that E]y and E^^ are not very correlated and that the bulk of the events is always in the 
small EiY^j) region, regardless of the value of E^{W). The peak in the double differential 
cross section is steepest in the jet scenario, however (see figure 30. c). 

Since Pt is conserved under Lorentz boosts along the beams, the W and the photon are 
exactly correlated in transverse momentum space in the jet scenario (i.e. Pt^/ = Pt^)^ so 
this correlation only makes sense in the 1 jet scenario, where it gets smeared and it peaks 
around the line Ptw = SOGeV regardless of the value of Ptj (see figure 31). 

The pseudorapidity correlation between W and 7 clearly shows the effect of the partonic 
zero as we can appreciate in figure 32. c in the jet scenario. The double differential cross 
section shows a valley in the negative 7]^ region along the whole range of 77^. Inclusion of 
1 jet processes (see figure 32. b) would push the valley up producing a plateau instead, as 
seen in figure 32. a for the 2 body inclusive scenario and smearing out the dip in the single 
distribution, as already noted in figure 11. 

Due to the exact correlation between Ptw and Pt-y in the jet scenario, the correlation 
between Ptw and r/^ in the jet scenario (figure 33. c) looks similar to the double differential 
distribution of Pt^ vs 1]^^ (figure 29. c), i.e., it shows an increasing effect of the partonic zero 
as Ptw decreases so that the photon rapidity dip becomes relatively deeper at small Ptw- 1 
jet processes present no dip in this double correlation, as seen in figure 33. b, and their effect 
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when looking at the 2 body inclusive process in figure 33. a is to smear out the dip, which 
reappears below Ptw = 30 GeV. 

The pseudorapidity of the W shows dips in all three scenarios (figure 6) and these dips 
are enhanced if Pt^ is constrained to be smaller than 30 GeV, as inferred from figures 34. a, 
b and c. We contrast this with the results for the rapidity of the photon, which only shows 
a clear dip in the jet scenario, even when looking at regions of low Ptj (figures 29. a, b and 
c). 

The Ptjet correlations with Pt^ and Ptw are quite different: in the Pt^ versus Ptj^t space 
the events are mainly concentrated along the small Pt^ and small Ptjet regions (see figure 
35) while in the Ptw versus Ptj^t space they are concentrated along the Ptw — Ptjet line 
(see figure 36). Something similar happens in the case of rjjet correlations with Pt^ and Ptw' 
in the Pt^ versus rjjet space the events are concentrated mainly along the small Pt^ regions 
(see figure 42) while in the Ptw versus rjjet space they peak along the Ptw = 30 GeV line 
(see figures 41 and 2). 

Although the single rapidity distribution of the photon does not present a dip in the 1 
jet scenario (figure 11), the dip in rj^ reappears if we fix the rapidity of the jet, as can be 
appreciated in figure 37. The dip moves up in 77^ as rjjet is increased and it is steeper in 
the regions of negative rjj^t- This figure also shows that the double differential cross section 
peaks along the line rj^ = rjjet- 

Figures 43 to 45 show the correlations between cone sizes Rij- Rjet,'y versus Rwjet (figure 

43) shows a symmetric shape around the plane Rwjet = tt at which the correlation peaks 
with a dip around Rjet,-) = 1; these features still show up in the single distributions in figure 
21. The correlation of Rwjet with Rw,^ (figure 45) peaks around Rwjet = Rw,'y = ti" and 
also has a symmetric shape with respect to the plane Rwjet = tt, confirming again the 
observations regarding the plot of Rwjet in figure 21. Although Rjet,'y versus Rwjet (figure 

44) shows no simple symmetry it also peaks near the planes Rjet,-y = and Rw,jet = ^r. In 
general the Rij cone sizes peak around Rij = tt and fall off quickly above Rij = 4, and, in 
particular, Rw,jet is symmetric around tt even when plotted against other variables. 
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V CONCLUDING REMARKS. 

We have completed an analysis of the exclusive reactions pp — > W'^'-f and pp — > W'^'j + jet 
by generalizing the methodology used by other authors in the context of heavy quark and 
Z pair production which consistently includes all divergent regions of phase space in 

the framework of the factorization theorem and dimensional regularization. We have taken 
into account order as QCD corrections and leading logarithm bremsstrahlung contributions 
and used the latest CTEQ set of parton distribution functions, which fit the newly released 
HERA data. We have treated the W boson as a real particle with a mass of 80.2 GeV. 

All our analytical results were presented in our paper [|I| and have been used in a Fortran 
code with Monte Carlo integration techniques. Histograms of single and double differen- 
tial cross sections and correlations for all outgoing particles {W boson, photon and, when 
applicable, jet) have been obtained in each of three experimental scenarios, namely 2 body 
inclusive, 1 jet and jet scenarios. 

Our results for total, single and double differential cross sections and correlations for the 
production of and photon accompanied by jet show the smallest theoretical uncertainty 
under variations of the mass factorization and renormalization scales when compared with the 
predictions in the 2 body inclusive and 1 jet scenarios. This means that accurate predictions 
are available in the jet scenario for quantities related to the W boson and the photon 
without including order a| corrections, however, they would have to be included for better 
accuracy in the other two scenarios. 

Previous work on the reaction pp — > W'^'-f + X has been devoted to the study of single 
photon distributions, photon- boson pair mass and photon-charged lepton pseudorapidity 
correlations. The complete set of distributions and correlations including the W boson and 
the jet therefore complements the studies of the Electroweak sector of the Standard Model 
(i.e. the magnetic moment of the W) and provides further checks for the QCD sector and 
the photon bremsstrahlung process. We note here that our results for the jet and 2 body 
inclusive reactions include leading logarithm photon fragmentation functions, so deviations of 
the observed experimental data from theoretical predictions could in part be accounted for by 
errors introduced by this approximation. To discriminate the effect of poorly known photon 
bremsstrahlung contributions from the effects of deviations from the Standard Model, the 
analysis of distributions and correlations in the 1 jet scenario, which, at the present order 
in perturbation theory are free of photon bremsstrahlung, would provide a valuable tool. 
With regard to total cross sections the two reactions pp — > W'^'-f and pp — > W'^'-f + jet are 
of similar importance, but, as we pointed out before, the latter reaction really requires even 
higher order QCD corrections to provide more accurate predictions. 
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